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Abstract. We investigate theoretically the Seebeck effect in materials close to 
a ferromagnetic quantum critical point to explain anomalous behaviour at low 
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by a spin- fluctuation temperature appears in the temperature dependence of correction 
terms as in the specific heat. 
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1. Introduction 

Experiments on clean materials near ferromagnetic quantum critical point (QCP) have 
revealed unusual properties, including non Fermi liquid transport and unconventional 
superconductivity. P, Wi The effects caused by quantum critical dynamics of spin 
fluctuations on the specific heat coefficient, the spin susceptibility, the resistivity, and 
so on, have been elucidated analytically at low temperatures. [SI IH |5l (Gj [7] In most of 
such theoretical analyses made so far, critical spin fluctuations are regarded to stay 
in thermal equilibrium. On the other hand, one may conceive of its inequilibrium 
counterpart of anomalous behaviours as well, which would be of fundamental interest 
too and should be paid due attention theoretically. As a representative of such 
phenomena, there are observations suggesting spin-fluctuation (or paramagnon) drag 
thermopower. In the Seebeck coefficient S{T) of UAI2, for example, there have 
remained a structure at low temperature, which is observed experimentally, |8l |9] but left 
unexplained theoretically. [T0| ITT] Among others, the most typical clear-cut experimental 
evidence would be those reported by Gratz et aL,[T2l[T3l[Tl] where the pronounced low- 
temperature minimum in S{T) of strong paramagnet RC02 (R=Sc, Y and Lu) was 
attributed to the paramagnon drag effect. Recently, Matsuoka et a/. [15] found a similar 
structure for AFe4Sbi2 (A=Ca, Sr and Ba). In effect, Takabatake et aL[T6] made it 
clear that the anomaly in S{T) is indeed caused by the ferromagnetic spin fluctuations 
prevalent in the materials by showing that those structure is completely suppressed by 
applying a uniform magnetic field. In contrast with the accumulating experimental 
evidence, there seems no theory to compare with the experiments available so far, 
with the exception of a brief account on a qualitative effect expected for localized spin 
fluctuations around impurity sites of alloys. [17] In this paper, we discuss an effect of 
uniform spin fluctuations in a translationally invariant system, and intend to provide a 
more solid footing on which to discuss the phenomenon. 

In section [2l we give an outline of a two-band model, which we adopt as a relevant 
model, along with approximations and assumptions conventionally made. In section [HI 
we introduce a function $f to represent inequilibrium displacement of spin fluctuations. 
In section IH we discuss that the leading effect of spin fluctuations appears on the 
T-linear term of S{T). In effect, in section 14.21 we discuss that the leading term 
contribution follows a universal relation to the specific heat, that is, q = eS/C ^ ±1 
revealed by Behnia et a/. [18] In the higher order terms, we have to consider not 
only a critical effect originating from equilibrium quantities, but also a genuinely non- 
equilibrium effect which has not been investigated before. In section [5], we investigate 
the latter contributions to find a characteristic temperature dependence, and the results 
are summarized in the last subsection 15.41 In section [6l we discuss the results and 
comparison is made with experiment. 



Spin-fluctuation drag thermopouier of nearly ferromagnetic metals 



3 



2. Model 

Let us introduce a two-band paramagnon model, which is conventionally employed 
to explain an enhanced resistivity of transition metals at low temperature. O |19l [20] 
The model has been applied successfully to explain, e.g., a saturation behaviour at 
elevated temperatures by taking into account a proper temperature dependence of spin 
susceptibility. [211 [22] 

The model is comprised of two types of electrons, i.e., wide-band conduction 
electrons and narrow-band itinerant electrons on the border of ferromagnetism. We 
denote the former as the s electron and the latter as the d electrons, representatively. 
The Hamiltonian consists of three parts, 

H = Hg + + Hd- 

The free Hamiltonian of the s electron is given by 

ka 

where c|,^ and Cka are the creation and annihilation operators for the electron with 
momentum k and spin cr. For simplicity, it is often assumed that the s electrons make 
a parabolic band with mass m^, i.e.. 

At each site i, they are scattered by the spin Si of the d electron at the same site through 
the Kondo s-d coupling, 

//,rf = J5]s,-5„ (2) 

i 

where J denotes a coupling constant, and Sj = ^J^aa' '^la'^aa'Cia is the spin of the s 
electron at the site i expressed in terms of the Pauli matrix vector To-o-' • Similarly, the d 
electron spin at the site i is given by 5, = | J^aa' dj^r^a'dia in terms of the creation and 
annihilation operators d}^ and c/jo- for the d electron. Spin dynamics of the d electrons 
is described by the Hubbard Hamiltonian, 

Hd = Y^ £d{k)dl^dka + UY^ rii^riii, (3) 

ka i 

where Uia- = d\^di„ [a =t, \) is the number operator of the d electron at the site i. 
The on-site repulsion U is fixed such that the d band is nearly ferromagnetic. To make 
analytical evaluation feasible, it is often assumed further that the d electrons are also 
parabolic with a different mass heavier than m^, i.e., 

and md 3> m^. The latter inequality is regarded as the basic ingredient of the model. 
Hence the d electrons act as heavy and fluctuating scatterers against the s electrons 
through the coupling of (|2]). In effect, this is taken into account as the second order 
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effect with respect to the couphng J, i.e., through the Born approximation. [T^ Then, 
the d electron comes into play through the (transverse) spin susceptibility In 
the random phase approximation, it is given by 

where 

jo _ jO 

Xo(9,c^) = Y1 I ^- (^) 

Y ^d{k + q)- ed{k) -uj -id 

Here, = /^(^^(A;)) = l/(e(^<<(^)~^)/'^ + 1) is the Fermi distribution function, and 5 is a 
positive infinitesimal. To investigate critical properties at low temperatures, [23] (E]) is 
expanded for small q and u/q as 

Xo(g, uj) = pfa - -^f + ^I^) ' 

for u < 2q, where q = q/kp^d and uj = uj/eF,d are the momentum and energy 
normalized by the Fermi momentum kF,d and the Fermi energy eF,d of the d electron. 
PF,d = ^dkF,d/'^T''^ is the density of states (DOS) at the Fermi level of the d electron per 
spin. Substituting (I7j) into ([5]), we obtain 

^^i' i-^^ ' 

for u < 2q, where tJ = PF,dU, and Kq = 1 — f7(^ 1) represents the distance to the 
QCP. 

The intrinsic transition probability Q^^'' that an s electron with momentum k is 
scattered to k + q hj absorbing a spin fluctuation with q and u via the coupling in ([2]) 
is given by 

Q 72 

2'rH = ^^(^'^)' (9) 

where S{q, u) denotes the Fourier transform of the spin density correlation function, 
which is related to the dynamical susceptibility by the fluctuation dissipation 
theorem. [23] 

2 

S{q,(^)= ^_ ^) ■ (10) 

The equilibrium transition rate is given by 

n't' = / dc.(l-r(e.(fc + g)))r(5.(fc))n°(a;)Q^+»5(a; + e.(A:)-e.(fc + g)), (11) 

where f^{es{k)) is the Fermi factor for the s electron, and n^(uj) = — 1) is the 

Bose function. With this g'', transport coefficients are derived by following the formal 
transport theory of Ziman[25 (cf. Appendix A). Transport properties of the s electrons 
in an electric field E and a gradient of temperature VT are described by the Boltzmann 
transport equation. 
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where Vs{k) = Vk£s{k) is the velocity of the s electron, and e(< 0) is the electronic 
charge. The right-hand side in f|T2|) is the collision integral for the s electron. 

To linearize the transport equation for the conduction electrons, a function $| to 
represent the displacement of the distribution function f{es{k)) from the equilibrium 
one f'^{es{k)) is introduced, i.e., by 

fie.ik)) = f%em - ^^^^n- (13) 

On the contrary, the d electrons are commonly assumed to stay in equilibrium, despite 
the applied fields. Then, for the collision integral in ffT2|) . we obtain 

q 

For definiteness, let the fields E and VT be in the direction parallel to a unit vector 
u. For the isotropic model, the magnitudes of the electric and heat currents due to the 
s electrons are given by 



L (15) 



and 



= 2 E ^ • vs{k){es{k) - /i) (-^^^^) n- (16) 

The factor 2 in front of the k sum accounts for the two spin components. As noted 
below (|T3l) . it is conventionally assumed that the corresponding currents due to the d 
electrons are neglected against the s electron currents. 

To obtain a solution $|, one may set ^f, = tu ■ Vs{k), while the constant r is fixed 
by the equation. Consequently, for the electric resistivity R = Ro(T) and the diffusion 
thermopower coefficient S = Sq(T), we obtain 



and 



where 



k.q 

The ordinary diffusion thermopower in f|T8|) is linear in T at low temperature, and is 
often expressed as 

QSfn.-, _ 7T^T d\ogas{eF,s) 
in terms of the spectral conductivity cTgi^) of the conduction electron. 
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3. Spin-fluctuation drag 

As remarked above, the d electrons are customarily assumed to stay in equilibrium 
regardless of the applied fields. To generalize the above framework to describe spin 
fluctuations with a shifted distribution theoretically, let us consider a bare dragged 
susceptibility Xoili^)) which is obtained by shifting uniformly the equilibrium bare 
susceptibility xoili^) in ([6]) by a small but finite amount Qq in momentum space. 
Similarly, we may define for the full susceptibility w) as well. Hence, 

X'^°{q,oj) is strongly peaked at q = go- 
First we derive a simple relation between xVili^) ^'^d Xfiili^)- According to 
([5]), we will obtain a similar relation for the full susceptibility. For the derivation, we 
introduce a shifted energy of the d electron, 

el'{k) = Sdik - go) ^ ed{k) - qo ■ v^ik), (20) 

where Vd{k) = Vk£d- Then, Xo^Q^^) is obtained by distributing the d electron with 
momentum k according to the shifted distribution that is to say, by 

Y ed{k + q)-ed{k)-u 

f\ ea{k)) - f^Edjk + q)) 
, \k + q)-e-/\k)-uj- 
Thus, by (120|) . we obtain the relation 

Xo°(?>^) - Xo(g,w + go ■ -^dlg))- (23) 

This is the result on which we base ourselves in the following. 

According to (|23|) . the drag effect is described by a function $^ = go ■ Vd{q). To 
understand what this represents, it is instructive to consider the isotropic case of (jl]), 
where Vd{q) = q/md- In this case, we obtain ^g = V-q where V = qo/md denotes a 
uniform drift velocity of the d electrons, or the spin fluctuations. In effect, the energy 
e~^'^°{k) represents the excitation energy of the d electron in the moving frame drifting 
with the velocity V . This is just a Galilean transformation. Indeed, noting that we can 
write 

d f{ed{k) )^ 
ded{k) 

and comparing this with (JT3i) . it would be clear that the new function $f represents 
the distribution shift of the d electrons, just as $| does for the s electrons. Thus, we 
argue that the drag effect of spin fluctuations is described in terms of in the way 
that Xdra.gil,'^) of the dragged fluctuations is represented as 

Xdrag(g,w) =X(g,W + $J), (24) 

in terms of the equilibrium susceptibility x(g,C(;). 

I A similar consideration was taken to derive the Drude weight of a Fermi Uquid. [25] 



/°(^^"(fc)) = rMA:))- "^.r..7^ < 
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Given the above argument, we have next to investigate how the formahsm in the 
last section should be affected by a no n- vanishing The first effect is to modify the 
collision integral in ( |T^ . To see this, here we follow how ( fT4l) is derived. The collision 
term in the right-hand side of (fT2|) is explicitly given by 

A = - E / dc^ [(1 - fk+g)fkn%u) - - fk){n\co) + 1)] (25) 

q •' 

y,Ql^^\u)5{u + es{k)-es{k + q)), (26) 

where we denoted fk = f{6s{k)) for the distribution function. According to the condition 
of detailed balance, the equilibrium distribution functions /° and n°(c<;) satisfy the 
relation 

(1 - - - fk)in» + 1) = 0. (27) 

Accordingly, by substituting ( fT3l) into (1261) . we obtain ( JT^ to the linear order in $|. To 
go further to take into account the inequilibrium shift of the d electrons, we regard that 
Qk^q'^iuj) in ( !26l) . or V^^'^ of ( fTTj) . depends on Xdrag(5', i^) in place of x{q,<^)- Then we 
can make use of ( l24l) . The first effect of $^ is to change the scattering probability "P^^'', 
which eventually has no effect owing to (1271) . The second is to replace n^{u) in (!26|) by 

n\uj - $2) ^ n\iu) - — (28) 
As a result, we obtain 

fk = -^i: {n + n^,) vz'- (29) 

q 

At this point, (!29l) clearly indicates a close analogy to the similar problem of phonon 
drag. [21] On the one hand, we can reproduce the previous results under the assumption 
$g = of no drag. On the other hand, owing to $g in fl29p . we can recover the correct 
identity fk = when the model is genuinely isotropic as implied by ([1]) and (jl]). In fact, 
in this case, we may set 

^l = ^i = u- fc, (30) 

where we put V = u without loss of generality. Then the null result for (129!) obtains 
from the total momentum conservation. This means that, if properly treated, the 
model should give no resistivity at all, irrespective of strong scatterings with spin 
fluctuations. In effect, the spin fluctuations in the inequilibrium state represented by 
(130|) are completely dragged along with the conduction electron currents. It is the 
fully dragged state in which all the s and d electrons drift with the same uniform 
velocity V, independently of the electric field E. This is the opposite limit to the 
case $g = without drag. In practice, in any case, we should have a finite rate fk 
by some mechanism neglected in the simple model, e.g., by Umklapp scatterings or 
by scatterings with extraneous agents. Moreover, generally, in order to investigate the 
degree of drag quantitatively, e.g., the temperature dependence through a wide range 
over a characteristic spin fluctuation temperature, $^ should be determined consistently 
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on the basis of its own transport equation. In general, the k dependence of $| and $^ 
may not be as simple as in fl30l) . 

For definiteness, however, we restrict ourselves to the low temperature regime, 
where we make use of the full drag assumption, ( !30!) . to elucidate non-trivial effects 
arising from our extra degree of freedom, A formal theory to discuss a general case 



is given in Appendix A 



4. Leading effect 
4.I. Limiting cases 

In the original model, the d electron currents are neglected on the basis of the basic 
inequality |fs(fc)| ^ \vd{k)\, or ^ vTis-IMUSH] Close inspection indicates that this is 
concluded through the additional implicit assumption $^ = u • Vi{k) {i = s,d) on the 
solutions of the transport equations, namely, by $| ^ $f ~ 0. As we saw above in 
(j30l) . this does not hold true in the presence of the d electron drag. In effect, the leading 
term contribution to the thermopower will arise from those dragged d electron currents, 
which would outweigh the normal diffusion term Sq{T) in f|T8l) due to the conduction 
electrons by a factor of ma/nis ^ 1. 

We obtain from ([1]), (fT5|) . and = u ■ k, 

2e 

j; = J,[$'] = Y^F<skF,spF,s, (31) 

where vf,s = kp^s/fris is the Fermi velocity. Similarly, ( !T6|) may be written as 

Ul ^ = TT^'ir^. (32) 

oe 0£ F s 

where Sp^s is the Fermi energy. The latter is obtained by expanding the integrand in 
( !T6|) with respect to the excitation energy es{k) — fi. The factor of vr^T^/S derives from 
the energy integral over es{k) to replace the k sum. Hence, from (ITS]) we obtain the 
ordinary T-linear Seebeck coefficient 

In the same manner, the d electron currents are evaluated. We may use 

J', ^ = 2eE^ ■ Mk) (- %X'(ir ) ^^^^ 

in place of ( !T5|) . and = Ud[^'^] as in f|T6|) . with which we obtain 

as in (133|) . Formally, this represents the diffusion thermopower due to the d electrons, 
as 5*0 does for the s electrons. Therefore, we should expect 

l^ol » l^ol, (36) 
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for Sq is proportional to the mass mj. Still, it is remarked that Sq in (1351) is not a 
directly observable quantity in general. In fact, from flA.14|) . the total thermopower is 
given by 

S, = (37) 

Therefore, on the one hand, in the conventional case without d electron drag, where 
l-^s I ^ l^dl |f/g I ^ \U^\, we recover the normal result 5*0 — S'g. On the other hand, 
in the opposite limiting case of the full drag, as the two currents and become 
comparable with each other, we expect a sizable modification from the normal result. 

To make this explicit, we remark that the currents are conveniently expressed in 
terms of their electron numbers and n^. In effect, it is straightforward to show 
Jg = 2nse from ( 13T]) . or more generally, we get it by a partial integration as follows. 

j:^2e^u...(,)(-'-g^)u.k (38) 
= -2ej;u. Vt(u ■*/"(£.(*))) +2e5; /"(£.(*•)). 

k k 

The first term represents the contribution from the Brillouin zone boundary of the k 
sum, which vanishes when the states there are unfilled. The second sum gives the result 
of the total number times e. Similarly, we obtain = 2nde for the d electron. These 
results simply represent that the whole electrons are drifting all together, as noted in 
the last section. Hence, from fl37|) we get 



Us + Ud 

Especially, in the limit rid ^ Ug, we obtain the enhanced diffusion thermopower Sq ^ Sq 
given in (l35l) . which is wholly due to the d electrons carrying the spin fluctuations. 

4-2. Equilibrium effect 

To the extent that we make use of an approximate expression ^n^e as above, 

one may obtain 11^ — 26^ correspondingly similarly, where generally represents free 
energy of the d electrons. Then we obtain 

S'o"~erf/(Tn,e). (40) 

This expression may be valuable as it is expressed in terms of the equilibrium quantities, 
which have been vigorously investigated. For example, one may have recourse to scaling 
argument for ed.[7\ We obtain S'g oc T by oc normally, while at the QCP, 
Sq oc TlogT"^ according to Cd oc T^logT"^ In terms of the electronic heat capacity 
C, one may substitute = CT to obtain 5*0 — C/{nde), or 

1^^ = -, (41) 

where n = + rid and S* ~ 5*0 — ridSQ/n under (l36l) . For hole like carriers, following 
as in fl38p . we find that the number rid becomes negative with the absolute value 
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\nd\ representing the hole number. Thus our drag mechanism supports the material- 
independent universality in q as revealed by Behnia et a/.|18] This is contrasted with 
the explanation by resorting to dominant impurity scatterings. [27] 

To go further to investigate the next order contributions, we have to consider not 
only those originating from the equilibrium quantities, which may be related to singular 
behaviour of the specific heat, but also the non- equilibrium effect which manifest itself 
in linear response to an applied field. The latter, though potentially important, has 
not been investigated before. In the next section, we focus ourselves to such singular 
contributions which vanish at zero field = 0. We find similar temperature dependences 
as that expected from the equilibrium effect through (HOll . 



5. Sub-leading corrections 

5.1. Extra currents 

The effect of spin fluctuations on the single particle excitation of conduction electron is 
described by a particle self-energy T,{k,e). The dragged spin fluctuations bring about 
a similar effect as those in equilibrium affect the thermodynamical properties. [3l H] We 
pay attention to the extra quasiparticle currents induced by the change of states at the 
Fermi level, as they are expected to make dominant contributions. We write an energy 
shift caused by a non- vanishing factor $f as 6es{k). Then the extra currents are given 
by 

J,[<l>^] = 2eE^ ■ ^^(^)^r(^'^^'^^^' ^^^^ 

and 

Us [$^] =2Y^u-Vs{k){e,{k)-fi) —L-6es{k) . (43) 

The effective energy of the conduction electron at the Fermi level is given in terms of 
the real part of the self-energy ReS(/i;,£:) by 

£^(fc) + ReE(M) 

^ s\ ) — 75 ■ 

1 - — ReE(A;,0) 
oe 

For the self-energy, we are interested in those part induced by the dragged spin 
fluctuations, which we denote as 5 (ReS(A;, 0)). Thus we have 

5es{k) ~ 5 (ReS(A;, 0)) + - (ReS(A;, 0)) , (44) 

as we need 5es{k) and 5(ReS(/c,0)) only to the linear order in The first and the 
second terms in (jH]) contribute mainly to Js[*!&'^] and f/^ respectively. In effect, we 
find 

U,[^^] =2Y^u- v,{k)i6,ik) - ^)2_^_<5 (ReS(A;, 0)) 
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-2{{v,{k)-u)—6{ReE{k,0)))k,, / p.(£)d£(£,,, - /i)^^ (45) 
oe ' Jo oEs^k 

PF,sI'm\ (46) 



r , , .2 df^ „ 

Jo oe.k 3 



3 

where /'(O) is the derivative at e = of 

I{e) = {u-v,{k)6{ReT.{k,e)))k^^,. (47) 

The angular bracket in fHTj) represents the average over the Fermi surface. In fHSl) . Ps{^) 
is the DOS per spin of the s electron, and pF,s = Ps{.^f,s)- Furthermore, we used 

''S,k 

Similarly as ( H6|) . we obtain 

JJ*'^] = -2p^,J(0), (48) 

using I{e) in f H7|) . As we find /(O) is insignificant, a correction to the thermopower due 
to the s electrons affected by the spin fluctuations is given by 

T( J| + Jj) 3e + rirf ^ ^ 

5.2. Self-energy 

We employ the self-energy in which a spin fluctuation excitation is emitted at one vertex 
and absorbed at the other one. It is given by 

S(fc, e„) = -\j'TY.Y. G{k', e'Jxik - k', - e'J (50) 

where e„ = {2n + l)7rT and e'^ = {2n' + l)7rT are the fermion Matsubara frequencies, 
G{k,en) is the temperature Green's function for the s electron, and x{Q,^n) is related 
to the d electron susceptibility ^(q,^) at the imaginary frequency u = iun, where 
Un = 2n7rT is the boson Matsubara frequency. By an analytic continuation, we obtain 
the following relation for the retarded functions, denoted below with the subscript R, 
which are analytic in the upper half plane of the complex frequencies; 

ReEnik, e) = —J^Y. I ^ImG/?(fc', w)ReXi?(fc -k',e~ u) tanh ^ 
2 ^, J~oo 271 2,1 

-Ij^Y. r ^ReG^(fc',^-c^)ImX/?(fc-fc',c^)coth^. (51) 

2 ^1 J-oo /TT 11 

To obtain the effect of we substitute x^(g, w) = Xdr^gili^ + i*^) from ( 124|) . Hence 
the shift 5 (ReS/j(fc, e)) is obtained from (ISTl) by substituting ^^^^[7^$^ in place of 
XR{q,(jj). For Gjiik'.uj), we use a free propagator Guik^uj) = 1/{uj — C,k + 16), where 
^k = £s{k) — /i. Owing to ImGij(fe', w) = —7i6{ijJ — C,k') and ([8]) for xr{q,^), the first 
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— ImGjj(A;', w)— RexR(fc - fc', £ - tanh — 

= E §-^^^XR{q, e - i,.,) (l - 2/0fe_,)) 

As this give only a convergent result, we neglect this part. Using fj30|) for $g, for fH7!) 
we find 



j2 



OO e — u + Vs[k) ■ q ou 2T 



2(27r)4 

where we substituted ^k-q ~ ^ i^s(^f,s) • 9, which holds in the important integral region 
of small \q\. Integrating over the angle between Vs{kps) and q, we obtain 



J2 .2kp^ 



2(27r)3 Jo 



q'^dq / dw 2 - 



e — u 



VF,sq 



log 



e — uj + Vp sQ 



duj 2T 



(53) 



In the parenthesis, only those terms odd in u contribute to the integral over u. Hence 
we find /(O) = 0, and the leading term in \u/{vF,sq)\ gives 



/'(O) 



Jo 



dq / cuduj- 



du 



1 

coth — - 
2T 



dq / wdw coth ■ 



2n^Vp^ Jo 



dco 2T 
To put in this expression, we may write the susceptibility in as 

— Imxij g,w = — - ^ . 

ouj du \q [k? + q^Y + {Coj/q) j 

for uj < 2q, where A = 36npF^d/U^ , C = Svr, and 

= 12KI/U = 12(1 - U)/U. 

We find 



(54) 
(55) 



/'(O) = (Xo +X(r)) = -36 f 4 



T\ 2 



where J = Pf,sJ, 



U \k 



(Xo+X(T)) 



and 
X(T) 



1 ^2kF,s/kF,d dq f"^^ _ _ d 
udu 



2 Jo 



UJ 



q Jo diu \{R^ + fy + {Cu/q) 



2kp^s/kF,d dq _ _ d 
u)du) 



q Jo dco \ (n^ + q2Y + ^Ccu/q) 



2 n [eF,d(^) 



(56) 



(57) 



(5J 



(59) 
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The former Xq is the part independent of temperature T, while the temperature 
dependence in the latter X(T) arises from the Bose factor n^{uj). In particular, for 
K = 0, we obtain 



Xq = - / d(7 / wdw 



{2kp^Jkp^af ( I 1 g4 

^^^7^7^ + ^ log ^ 



^ {kFjkF,dY + C^- 
5.3. Temperature dependence: X(T) 

To obtain an explicit expression for the temperature dependent part 1.{T), we adopt an 
approximation to set 

n (u) = < ^ 61 

^ ' \ 0, u > cT ^ ' 

where c is a constant of order unity (cf. ( IB. 21) ). Consequently, we obtain 

X(T)=X„(T)+Xb(T), (62) 

where 
and 

1 / T y r{2kF,s/kF,if ^2 f 9 ( u \ 

^'^^^ ^ 2 \7^d) kl ^'^ Jo "^""d^ («^' + 9')' + V ■ 

Here we introduced a characteristic scale for the normalized momentum, 

cT , , 

^0 = (65) 

We may take the limit k = for fl63p to obtain 

rm 2T ^go dg gpT T ^ 



^2 



On the other hand, for fl64j) . we obtain 



c ( T Y r(2kF,s/kF,d? dg 



for which the main contribution comes from around the lower limit of the integral. Let 
us discuss two cases depending on the relative size of R and separately. 
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First we consider the case H/qo ^ 1, which is the low temperature hmit for R > 0. 
In this case, we obtain 

c f T Y 1 /•(2fcF,.AF,d)' df 



^c/n^i,„,(i^^Zin^)!±(c*Z5!)!. (68) 

In terms of a characteristic temperature of spin fluctuations defined by 

C TlU ' ^ ^ 

we find 

UT) ^ ^ Of log / ^ ■ (70) 

In the literature, a spin fluctuation temperature, 

Tsf = £^,rfiro' = £i.,,(l-[/), (71) 

is commonly used as well. Indeed we have T ~ Tgf for ~ 1. Lastly, in the quantum 
critical limit /t/go ^ 1, we obtain 

2 



1 / T \ /•(2fcF,s/fcF,d)2 dg2 



9(2Cgo)4/3 1,£^,J ~ 3y3(Cc)4/3 1,£^,J 



(72) 



5.4- Results 

We may neglect TaiT) against T^iT), for T <^ 6^7^^. For fH9l) . we obtain 

A^, ~ ASo^ + A^,(r), (73) 

where 

° 3e(n, + n,) l^f/yl 1^^^,^ ^ ( W^f.^)' + (Svr)^ ' ^ ^ 

and 

AS,(r) = ^i^f|)V|^)*P,,i,(r)r. (75) 

e[ns + nd) \U J \kF,s J 
The former ASq to modify the T linear term may be effectively neglected, while the 
latter ASs{T) gives a sub-leading correction. At the low temperature T <t^T, with ( 170|l . 

we get 

" skt;^ it/j ii;;:] "-^ (?) ^;^)^ • (™' 
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where we set c ~ 1 for simplicity (instead of (lB.2p ). In the opposite hmit, from fl72p . we 
obtain 

In the same manner as f/^f^*^] discussed above, one can think of an additional heat 
current Af/(i[<l>'^] caused by the intraband many-body effect due to the on-site repulsion 
U in the d band. Formally, the corresponding results are obtained straightforwardly by 
replacing kp^s, Pf,s, and 3J^/2 in the above results by kp^^, PF,d, and U'^, respectively, 
i.e., 



4 /T\2 4+(t/tY 

^^'^(^) ^ Q f ^ , PF4T 7f log , ^ . / , 
9e{ns + nd) \T J (t/t) 



(78) 

in place of fl76|) . The results are modified in some ways in generalizing the model. The 
constant of 4 in the logarithm of f l75]) stems from {2kp^d/ kF,d)^, where 2kF,d sets the upper 
cutoff for the momentum q of spin fluctuations. If we should have set a cutoff parameter 
Qc differently, the factor should be replaced by g^, where (jc = qc/kF,d- Moreover, if we 
had assumed a phenomenological coupling g between electrons and spin fluctuations 
instead of U, the results will be reduced by a factor of {g/U)"^. 

6. Discussion: comparison with experiment 

To compare the theoretical result S{T) with experiment, some assumptions like the free 
dispersions in ([1]) and @ should not be taken literally. In particular, the T-linear terms 
Sq {i = s,d) for ^0 would be able to have either positive or negative sign, depending 
on the energy dependence of the respective DOS at the Fermi level, while the relation 
|S'g I ^ \Sq\ will always hold true for their relative magnitudes. Therefore, as the leading 
effect at low temperature, we generally expect an enhanced T-linear term. 

So ^ S'o = ^^^0, (79) 

unless Us ^ n^. Effectively, this term is indistinguishable from the diffusion term 
contribution, as discussed below ( l35l) . It is indeed due to the drag current of the heavy 
d electrons. Without drag, we recover the conventional result ~ S'o — of the 
diffusion thermopower due to the conduction electrons. We expect that the latter holds 
true at high temperature T > T where the s-d scatterings become too weak to sustain 
the d electron drag. Therefore, it is reasonably expected that we should find some 
structure in the temperature dependence of the thermopower S(T) around T ^T, which 
is brought about by the crossover between the T-linear terms with different magnitudes 
of coefficients. This is schematically shown in figure [H 

Takabatake et al.^6\ have shown experimentally by applying the magnetic field of 
15T that an S-shaped structure in S(T) of CaFe4Sbi2 observed at low T < T 50K is 
suppressed to yield a normal T-linear diffusion term. This is consistent with our result 
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S 




Figure 1. The effect of spin-fluctuation drag on tlie tfiermopower S{T) is 
schematically shown. The bold lines are drawn to interpolate the two linear relations, 
S{T) — Sq and Sq. The low temperature Sq in ([75]) is due to drag of those heavy 
electrons pertaining to the spin fluctuations, while the normal diffusion term 5 ~ 
due to light conduction electrons appear at high temperature T T ~ Tgf, a spin- 
fluctuation temperature. The left (a) is for Sq > 0, while (b) for Sq < 0. The latter is 
compared with RC02 (R=Sc, Y and Lu) by Gratz et a/ . [HI [HI [Tl] 

for Sq < 0, Sq > and |5'q/S'q| ~ 8. In this case, the conduction band for Sq consists 
mainly of 5p states of antimony. Moreover, they have shown that the temperature 
dependence of the spin-fluctuation contribution AS" = S — Sq is not monotonic. To 
explain this theoretically goes beyond the scope of this paper, as it requires us to solve 
the transport equations concretely. Similarly known before were the low temperature 
minima in the thermopower of RC02 (R=Sc, Y and Lu), which had been stressed by 
Gratz et a/. [121 [131 E] cis the experimental evidence of paramagnon drag. Their results 
can be compared with our result for < 5^ < in figure [11(b). 

On the correction terms, it is generally expected that the d electron contribution 
ASd will become more important than ASs when the d electron current becomes relevant 
indeed. As discussed in section 14. 2[ we have to consider two sources of contributions, 
one due to the equilibrium effect and the other due to the non-equilibrium effect in 
section [5l Interestingly, we find that both give the same temperature dependence, 
S{T) = aT + /3T^logT^^ away from the QCP. Nevertheless, we notice an important 
difference. While we observe /3 oc 1/Kq from the results of the last section, /3 expected 
from a correction term in (l40l) has an extra factor of 1/-K'q.[3J This means that the 
equilibrium effect becomes more important. We suspect that this would hold true at 
the QCP too, although there has been no definite calculation deriving the corresponding 
free energy correction oc T^/^ in accordance with our result. 

In any case, we remark that the relative magnitude of the electron numbers 
and rid may have an effect on the correction terms, the sign of which will depend on 
the factor en = e{ns + Ud), that is, the direction of the net current. In most cases 
where the model applies, the current carrier in the heavy-electron band will be hole 
like. Moreover, we generally expect that \ns\ will not exceed \nd\, or the net current 
would be hole-like, e{ns + rid) > 0. Accordingly, ASi > (cf. (HTl)). This is consistent 
with a model calculation of the spin fluctuation effect on the resistivity, where Jullien 
et a/. [21], [Tl] pointed out the important role of the parameter ^ = kp^c/kp^ on the 
transport properties of spin fluctuations systems. We observe the dependence in our 
results of ( 1761) and ( 1771) . To compare their numerical results with experiments, they 
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10 20 30 40 50 60 
T[K] 

Figure 2. The points are the experimental data of S{T) for AFc4Sbi2 (A = 
Ba, Sr, Ca).[15^ The lines are the least squares fits by the theoretical expression in 
(|80l) with the parameters given in table [TJ 



should choose ^ < 1 generically, that is, \ns\ ^ \nd\. 

To conclude, let us fit the low-temperature experimental data for S(T) of AFe4Sbi2 
(A = Ba, Sr, Ca) reported by Matsuoka et aL[T5] with 

S{T) =aT + (3TU) log . ^ . / , (80) 

where a, /3, T and 6 are regarded as parameters. In table [H we present the fitting 
parameters obtained for 5 = 4 by the least squares fits of the low temperature part 
of the data for T ^ 17K (< T). The results are shown in figure |2l along with the 
experimental data points. We find that a's do not depend much on the other parameters, 
and the ratios of /3 and T between materials are nearly independent of 6. The relatively 
large values for a and /3 will be more properly ascribed to the heavier d band than to 
the conduction band, in accordance with our result. Note that these coefficients are 
susceptible to the equilibrium effect of mass enhancement, |3l H] which we did not take 
into account explicitly (cf. section l^^ l§l The positive /S implies that the net current 
is in the hole-like direction. The relative material dependence of (3 in table [1] may be 
qualitatively compared with the observed static uniform susceptibility xo pd/K^, that 
is, Xo(BaFe4Sbi2) : Xo(SrFe4Sbi2) : Xo(CaFe4Sbi2) ~ 1 : 1.6 : 2.5. 
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§ Owing to prevalent anharmonic phonons in this skutterudite system, it may not be a simple matter 
to extract the electronic contribution from the observed specific heat coefficient 7, which does not 
depend sensitively on the divalent ion A.|15) 
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Table 1. Parameters to fit the tlicrmopower of AFe4Sri2 in figure [21 







/3 [^,Y/K^] 


T [K] 


BaFe4Sbi2 


-0.35 


1.5 


48 


SrFe4Sbi2 


-0.74 


1.6 


37 


CaFe4Sbi2 


-1.4 


3.7 


58 



Appendix A. Formal transport theory 

The formal expressions for resistivity and thermopower referred to in the main text are 
derived by adapting a general variational method of Ziman,[23] according to which $^ 
in f lT^ and in ([23]) are regarded as variational trial functions. Below we substitute 
r7j$* for $* {i = s,d), and take variation with respect to the arbitrary parameters r]i. 

On the one hand, the microscopic entropy production rate corresponding to (l29ll is 
given by 

k,q 

= f H P^m■ (A.i) 

i,j=s,d 

The components of the matrix Pij defined in (lA.ll) are explicitly given by 

k,q 

Ps, = p,s = ^T. {n - n+q) ^t, (A.3) 



'P k,q 

Pm =^E<'K)'- (A-4) 

k,q 

In ( lA.ip . not only emission of a paramagnon corresponding to ( l29ll . but the reverse 
absorption process is also taken into account. In the special case of the full drag without 
Umklapp processes, there holds the relation — $| = $q by ( l30l) . so that we get the 
following identities, 

Pss = Pdd = —Psd- (A-5) 
On the other hand, the macroscopic entropy production is given by 

^macro = + U-V^. (A.6) 

In the linear response regime, the electric current J and the heat current U are written 

as 

J =r],J[^'] + r]dJ[^% (A.7) 

U = r^sU[^']+7^dU[^% (A.8) 

where J[$*] denotes the current flow caused by {i = s,d), i.e., J[$*] formally 
represents the part of the total current which depends linearly on [/[$*] is similarly 
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defined. In general, these currents have different functional forms. It is remarked that 
is not to be identified with the current in the i band. Owing to the interband 
interaction, the distribution shift in the i band can induce a current in the other 
band. 

The variational parameters rji are determined so as to maximize S'scatt after equating 
S'scatt and 5'macro-[2^ Substituting the solutions into (IA.7|) and flA.8|) . we obtain 

j= y: jm{p-%[jm-E-^u[<^^]-vTy (a.9) 



and 



U= J2 Um{p-%(j[<!>^]-E-^U[^^]-VTY (A.IO) 

i,j=s,d 

where {P~^)ij is the inverse matrix of Pij. For definiteness, let the applied field E and 
VT be in the direction of a unit vector u. In an isotropic system, or in cubic symmetry, 
the results are expressed with the magnitudes J[$*] = J[$*] ■ u and = ■ u. 

From ( ]A.9p . we obtain the electrical conductivity, 

l,m=s,d 

The resistivity R = is given by 

PsdPds 

R = Ro{T) -^^^ , (A.ll) 




where 

MT) = jj^. (A.12) 

The latter, given in f|T7|) . is the resistivity that we obtain when we have no spin- 
fluctuation drag. In fact, this is the central formula to explain an enhanced resistivity 
of a spin fluctuation system due to normal scattering processes with long-lived spin 
fluctuations. [T9 | [20 | [5] According to (lA.lip . the d electron drag modifies the resistivity 
in two ways. First, we note that the numerator in (lA.lip vanishes in the full drag case, 
flA.Sp . This represents physically that a finite resistivity is brought about only with 
those scattering processes which can degrade the total net current. On the basis of a 
more realistic model, a proper treatment of Umklapp scattering processes could make 
the numerator a non- vanishing factor of order unity. Secondly, the positive factor in the 
denominator has the effect of suppressing the resistivity. This is due to an additional 
drag current of the d electrons. When fully dragged, the d electrons carry n^/ns times 
as large current as the s electrons, where Ud/us is the ratio of the electron densities. In 
general, this would not be negligible quantitatively, and it might be so even qualitatively. 
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From the condition of no heat flow U = for (lA.lOp . we obtain the Seebeck 
coefficient, 

j[$'](p-i);™f/[$"^] 

S=^b!I^ . (A.13) 

l,m=s,d 

From this we can obtain the result for the the full drag case of (lA.SP formally as a 
special limit. It is expressed simply by the ratio of the total energy current to the total 
momentum current as 

_ 1 um + [/[$^] 

Indeed, the simple result in this limit is straightforwardly generalized to many-band 
models. It is owing to this simple property that we investigated this limit devotedly in 
the main text. 

On the other side, the case without drag is obtained for [/f^*^] = Jf^*^] = as 

- fA 15) 



as presented in (fTSj) . As a matter of fact, the no-drag results of f lA.12|) and ( ]A.15|) are 
directly derived without taking into account from the beginning. 

Appendix B. Temperature dependence of X(T) at k = 

To evaluate X(T) in ( l59l) . we made the approximation as given in ( l6Ti) . We obtained 
(!72|) for Th{T) in which signifies the main correction term of AS* oc T^/^ in the 
quantum critical regime. The exponent 5/3 is the same as for the resistivity. [5] The 
derivation in section 15.31 indicates that the important contributions come from a; ~ T. 
In effect, this is the upper limit of the u integral for q <^ %, and the high-energy cutoff 
is naturally provided by the Bose factor nP{u) in the integrand, without employing the 
approximation in (!6T|) . With this in mind, we can obtain the result for R = directly by 
transforming the integral and taking the limits for the bounds of integration as follows. 

I,(T) = ^ / df udu^ 2 ^ (^) 



2^2/3 i.2kp^Jkp^i/Ty-^ foo d 

dv 



u 



2e\^ Jql/TV^ Jo + {Cu/6F,d) du - 1 

2^2/3 ^oo ^oo d / U 

dv 



2e\^ Jo Jo + {Cu/eF,df du Ve« 

/ \ 2/3 , 8 / V 



3\/3C4/3 \ef4 Jo du 
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= 1.10— ^(^Y\ (B.l) 

By comparing fIB.ip and fl72l) . we obtain c~'^^^ ~ 1.10, or 

c ~ 0.928. (B.2) 
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